AD-A042  376  CALIFORNIA  UNIV  SAN  DIEGO  LA  JOLLA  DEPT  OF  PHYSICS  F/G  20/12 

OSCILLATORY  INDIRECT  INTERACTION  BETWEEN  ADSORBED  ATOMS. (U) 

JUL  77  K H LAU*  W KOHN  N00014-76-C-0050 

UNCLASSIFIED  TR-40  NL 


I OF  I 
A0AO4237e 


AD  A 0 4 2 3 7 6 


f 


TECHNICAL  REPORT  ^ 40 


CONTRACT  NONR-N00014-76-C0050 


PROJECT  NR  017-653 


Supervisor:  Professor  Walter  Kohn 


DEPARTMENT  OF  PHYSICS 
UNIVERSITY  OF  CALIFORNIA  SAN  DIEGO 

La  Jolla,  California 


P O o 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the  United,  \ 
rt  States  Government . ^ 


DISTRIBUTION  OF  THIS  DOCUMENT  IS  UNLIMITED 


July  1977 


• 


'r> . > 


1 


* 

Oscillatory  Indirect  Interaction  Between  Adsorbed  Atoms 

K.  H.  Lau^  and  W.  Kohn 
Department  of  Physics 

University  of  California,  San  Diego,  La  Jolla,  California  92093 


Abstract 

The  oscillatory  part  of  the  indirect  interaction  between  two 

adatoms  on  a metal  surface,  separated  by  a large  distance  R,  is 

examined  for  the  jellium  as  well  as  for  the  tight-binding  model.  The 

asymptotic  form  of  the  interaction  is  derived  from  the  singularity  of 

its  Fourier  transform  at  wave  number  q = 2 . For  a substrate 

with  a spherical  Fermi  surface,  the  interaction  energy  is  proportional 

to  cos(2  k R)/R  . The  fall-off  is  faster  than  the  corresponding  oscil- 
F 

latory  interaction  between  two  impurities  in  the  bulk,  which  varies  as 

'I 

cos(2  k R)/R  . For  substrates  with  a partially  filled  surface  band, 

F 

the  interaction  energy  falls  off  much  more  slowly.  Under  simple 
assumptions,  the  interaction  behaves  as  cos(2  k^R)/R  . Compari- 
sons are  made  with  the  numerical  calculations  by  Einstein  and 
Schrieffer,  The  experimental  results  of  Tsong  on  Re  atoms  adsorbed 
on  W using  field-ion  microscopy  is  also  discussed. 

Supported  in  part  by  the  National  Science  Foundation  and  the  Office 
of  Naval  Research. 

^lUM  gra<luate  fellow. 


i 

i 


i 


li 


1 


\ 


1.  Introduction 


In  earlier  papers,  we  examined  the  substrate -modified  dipole - 

dipole  interaction^^  ^ and  the  phonon-mediated  elastic  inte raction ^ 

(3) 

due  to  two  atoms,  A and  B,  adsorbed  on  the  surface  of  a metal, 

and  separated  by  a large  distance  R.  Both  of  these  are  non-oscilla- 
- 3 

tory  and  vary  as  R 

The  oscillatory  indirect  interaction  between  adsorbed  atoms 
has  been  discussed  earlier  by  Grimley,  Grimley  and  Walker,  ^ 
and  Einstein  and  Schrieffer,  and  Tsong^^^  reported  experimental 
evidence  for  the  oscillatory  interaction  between  Re  atoms  adsorbed 
on  W.  These  papers  will  be  discussed  in  Sec.  V. 

An  exact  expression  for  the  static  density  response  function 

/Q  \ 

K(r,  r')  for  a sine-wave  gas  has  been  derived  by  Rudnick.  (The 
electrons  are  confined  to  the  half-space  z s 0 by  an  infinite  potential.  ) 
The  response  function  is  given  by 


k(£.  e')  - 


\ K R.  o 1 ' 


;i.  1 ) 


where  k.^  is  the  Fermi  wave  number. 
F 


( 4 Me  ^ k . 


^TF  ~ / T~  Thomas-Fermi  wave  vector, 

nti 


^o  " -f  (y  - y')^  + (z  - z ' 1 ^ 


Z 


R = [(x  - x'  + (y  - y ’ + (z  + z ' ^ 


and  j^(x) 


sin  X cos  X 


is  a spherical  Bessel  function  of  first  order. 

We  have  expanded  this  function  K(r,  r^  for  the  case  where 


both  r and  r^  are  near  the  surface  and  the  separation 


R = f(x-x^)^  + (y-yM^]^^^  is  large  compared  to  k ^ . After  some 

r 


straightforward  algebra,  one  finds  the  following  expansion  in  inverse 
powers  of  R: 


z“^z  ' ^ cos(2k  R) 

£' ) = i " ° 

R \R 


(^) 


(1.  ;m 
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This  is  in  contrast  to  the  R dependence  of  the  Ruderman- 


Kittel  type  interaction  between  two  point  imperfections  in  the  bulk  of 

(9) 


a free  electron  gas. 

In  this  paper,  we  also  examine  the  indirect  interaction  between 
two  adatoms  under  more  general  assumptions.  In  Sec.  II,  we  use  the 
jellium  model  for  the  metal  surface  in  which  the  electrons  extend  be- 
yond the  jellium  edge,  and  in  Sec.  HI,  we  use  the  tight-binding  model 
for  both  the  metal  and  the  adatoms.  We  shall  show  that  the  P'ourier 
transform  of  the  interaction  energy  has  a singularity  at  wave  vector 


|q  I = 2 k , where  q is  parallel  to  the  substrate  surface  and  2 k 


, ^ --  ^ _ - ^ 

is  the  appropriate  diameter  of  the  Fermi  surface,  and  that  this 


'I 


w 


singularity  is  responsible  for  the  oscillatory  interaction  at  large  dis- 


tances. If  we  assume  that  the  metal  has  a spherical  Fermi  surface, 

then  we  find  that  the  interaction  energy  is  proportional  to  cos(2k  R)/R^. 

F 

In  the  vicinity  of  a surface,  the  amplitude  of  the  relevant  wave -func- 
tions tends  to  zero  as  R tends  to  infinity.  This  weakens  the  singu- 
larity of  the  response  function  compared  to  that  for  the  bulk  electron 
gas,  and  causes  the  interaction  to  fall  off  more  rapidly. 

In  Sec.  IV,  we  consider  the  case  where  the  substrate  has  a 
partially  filled  surface  band.  Virtual  excitations  and  de -excitations 
of  the  substrate  electrons  in  the  surface  band  contribute  to  the  inter- 
action energy  which  falls  off  much  more  slowly,  so  that  it  will  be  the 
dominant  part  of  the  interaction  between  the  adatoms  when  they  are 

far  apart.  If  we  assume  that  the  Fermi  surface  for  the  energy  band  is 

-2 

cylindrical,  then  asymptotically,  the  interaction  falls  off  as  R 

II.  Interaction  Energy:  Free  Electron  Model 

We  consider  a non-interacting  electron  gas  in  a jellium  back- 
ground, occupying  the  space  z ^ 0.  The  adatoms  A and  B,  at  r 

A. 

and  Tg  respectively,  are  represented  by  an  external  potential  of  the 
form 

ACP(r)  = . ( 1 r - r I ) + cp  ( I r - r I ) . (2.1) 

The  correction  to  the  ground  state  energy  of  the  electron  gas  is,  by 
second -order  perturbation  theory. 


4 


1 A"-P  I 0 ) 

E - E. 
o 1 


(2.  2) 


The  mutual  interaction  between  atoms  A and  B,  AE.  . , is  obtained 

int 

by  subtracting  from  (2.  2)  the  corresponding  expression  in  which  A 
and  B are  infinitely  far  apart: 


^ ( 0 |cp  g 1 i ) ( i 1 I 0 ) 


E - E. 
o 1 


+ c.  c. 


(2.  3) 


where  |0)  is  the  ground  state  and  |i)  the  excited  state  of  the  elec- 
tron gas, 

(a)  Electron  wavefunction 

The  normalized  ground  state  |0)  and  the  excited  states  (i) 
are  Slater  determinants  of  one-electron  wavefunctions 


(r)  = (4n^)  e^S  S k ). 

K Z 


(2.  4) 


p and  R are  two-dimensional  vectors 


p = (k  , k ) , R (x,  y ) , 
^ X y 


(2.  5 ) 


and  f(z;  k ) satisfies  the  equation  (in  atomic  units). 


^ dz^ 


f(z;  k ) = 0 , 
z 


(2.  6; 


where  v ,,(z. ) is  the  effective  potential, 
eff 


V ..(z)  vamshes  in  the 
eff 


interior  of  the  metal  and  rises  rapidly  to  a constant  value  A0  outside. 


j 


is  a point  near  the  metal  surface  such  that,  effectively,  v = 0 

for  z ^ z (Fie.  1 ). 
o 

We  shall  consider  only  small  values  of  k , since  these  give 
the  main  contribution  to  (2.  3).  In  the  interior  of  the  metal,  we  have 


f(z;  k ) - sin(k  z - y (k  ))  , 
z z z 


(2.  7) 


where  y{k  ) is  the  phase  shift.  For  z > z , we  can  drop  the  term 
■z.  o ^ 

1 2 

— k from  (2.  6)  since  it  is  small  compared  to  A0  , and  obtain 


f(z;  k ) = A(k  ) f(z)  . 
z z 


(2.  8) 


A(k^)  is  chosen  so  as  to  match  both  the  phase  and  the  amplitude  of  the 
interior  function  (2.  7)  at  z = z and  f(z)  is  the  solution  to  the  equa- 
tion 


(2.  9) 


with  the  boundary  conditions  f(z)  =-  0 as  z ^ and  f(z)  - z + constant 

for  z near  z . With  this  choice  of  f(z),  we  find  that  A(k  ) is  equal 
o z 

to  k , so  that,  for  small  k and  z > z 
z z o 


T I \ It  3.- 1/2  ip  • R 

(r)  = (4tt  ) e V A/  k f(z)  . 

k A'  z 


(2.  10) 


The  matrix  element  (i  | cp^  | 0)  , where  a = A,  B,  |0)  is  the 
ground  stalt;  of  the  system  and  (i/  the  excited  state  consisting  of  an 
elec  tron-hole  pair  of  wave  vectors  k - (p,  k ) and  k^  = (p^ , k^  ) is 

^ aj  Z ^ /V  Z 

given  by 


6 


1 


<i|cpJO>=-^f 
4n  J 


d^r  ■ S f(z;k' )cp  (VIR-R  \^+(z-z 

— ' /X/  ^ 


z a ' ~ ~ a' 


~l’(z;k  ) 
z 


1 iq  • R 

1 ^ ~a 


4tt' 


/dzf(z+z  ;k  )f(z+z  k' ) I d ReS  l/l 

J a z a z J a 


where 


, 2 Z, 
R +z  ), 


(2. 11) 


q = p - p 

/V  /^^  /V 


(2.  12) 


The  energy  denominator  in  eq.  (2.  3) 


E -E.  = --^(q^+2p*q  + k'^-k^ 
oi  2 Z ^ z z 


1 2 

vanishes  for  p - q = - — q and  k = k'  = 0,  so  that  the  main  contribu- 
Z Z c z z 

tion  to  (2.  3)  comes  from  p,  k and  k'  near  these  values.  Hence  we 

/V  Z Z 

shall  consider  only  small  values  of  k and  k^  and  write 


< i |cp  I 0)  = — ^ e^S'  ■S  k k'  1 (z  ) , 
' a ' ^3  '7.  7.  n n 

4rr 


z z a a 


(2.  13) 


where 


CX3 

/■ 


I (z  ) = / dz(f(z  -fz  ))  / d R e ^ ^ c;:  (vR 

a a / a I o. 


'/■ 


..iq  - R.  ,7p^  ^ 


(2.  14) 


The  inutual  interactit>n  between  adatoms  A and  B is  then 


r 2 ‘ 

q~^)  'A<“A"B"B’y‘'  3 “ ~ ” G(q)  + c.c. 

(2.  15) 


i-^) 


J 


where  the  factor  Z takes  spin  int(j  account  and 


G(q)  = - 


J f 

2.2  J 


dk' 


k^k'^ 
z z 


z F 

k ^ 0 
z 


^ -^(q^^+Zk  q+k'^-k^) 
2 X z z 


(2.  16) 


A cutoff  k has  been  introduced  in  the  integral  to  avoid  divergence, 

but  as  we  shall  see  later,  the  precise  value  of  k is  immaterial. 

^ o 

(b)  Evaluation  of  the  integral 

The  asymptotic  form  of  the  interaction  is  derived  from  the 

singularity  of  G(q)  at  2 k . Thus  we  shall  retain  only  the  part  of 

F 

G(q)  that  is  non-analytic  at  q = 2 k and  drop  the  remaining  part. 

F 

The  singularity  of  G(q)  is  very  weak.  In  fact,  G(q)  and  the  first 

three  of  its  derivatives  are  continuous  at  q = 2 k^  and  only  the 

F 

fourth  and  higher  derivatives  are  singular.  In  order  to  extract  this 
singularity  more  conveniently,  we  differentiate  G(q)  with  respect  to 
q four  times. 

For  q>  2 k , the  denominator  is  always  positive.  Differ- 
F 

entiating  G(q)  four  times  and  integrating  over  dk^  and  dk^  we 
obtain 


■G(q)  = 


dq 


J ^ " 


8(q+k  ) 

X 


6 (q  + k ) 

X 


2 .225  2 /224 

(q  +2k  q+k  -k  ) (q  +2k  q + k -k  ) 

X z z X z z 


1 


(q^+2k  q + k'^-k^)^ 
X z z 


Idk  dk  dk  dk' 
X y z z 


8 


= -6n 


^2  ,1/2 

- ,c,^2Kq-k^’'^'„^2kq-kV'^ 

L ‘ X z X z 


— ; 7T77 

/ 2 z 3/Z  X z 

(q  + 2k  q - k ) 

X z 


(2. 17) 


The  cutoff  k is  no  longer  necessary  since  the  integral  is  now  con- 
o 

vergent. 

For  q « 2 k the  main  contribution  to  (2.  17)  comes  from 
r 

k -k  and  k 0.  Putting  k = P cos  0 and  k = P sin  ^ and 
X F z X z 

expanding  around  9 = tt  , the  first  term  of  (2.  17)  becomes 


- 30tt 


^F  ^ 


4 3 

d0  (q  + P cos  9 ) P sin  9 


(k^  - 

F 

Z Z Z 5 / Z 

(q  + 2Pq  cos  9 - P sin  9 ) 


/42  21/23/  u du 

dP(q-P)^(k  -P  ) ' P^  / r . . . 2, 

^ / (q^-2Pq+(qP-P^u^f 


'I 


■-.I 


sin  P d 0 


J (q  - Zk  cos  0 ) 
o F 


,4rTT  , ,-1/2  n 1 

Nkj,  8kj.  4kj, 


(2.  18) 


for  q ~ . Hence,  for  q =»  2 k^,  , the  singular  term  in  (2.  18)  is 

1 2,5/2  -1/2 

~ kp  (q  - 2k^,)  . The  other  two  terms  in  (2.  17)  give  terms 

which  are  less  singular  than  the  first  and  can  be  neglected. 

For  q < 2 k^^,  , we  first  take  the  principal  value  of  the  integral 

over  dk'  in  (2.  16)  to  obtain 
z 


I z o 


dk  dk  dk 
X y z 


/ 'y  'y  \ I 'y 

— / k (q  + 2k  q - k ) dk  dk  dk  . (2.  19 ) 

2 J ? xz  xyz 

q'^+2k  q - k s 0 

Then,  differentiating  G(q)  with  resp<;ct  to  q and  proceeding  as  be- 
fore, we  find  that  the  leading  term  in  (q  - 2k  ) is  also  less  singular 


Thus,  the  most  singular  te;rm  of  G(q)  is  given  by 


_^G(q)  = - |q  - 2k^| 

dq 


R(q  - 2kp)  , (2.  20) 


X < 0 


X > 0 


(2. 21) 


For  large  |R  - R_|,  we  can  evaluate  the  integral 
^ A ^ B ” 


L 


where 


f , iq  • (Ra-Rr)  o •,  cos  Zk  R) 

/d  q e G(q)  ^ - n k i O 

J - ^ VR 


ii) 


(2.  22) 


R = ) R . - R_1  , so  that  the  interaction  energy  is  given  by 


AE 


cos(2  k R) 
IJz  J I (z„)  ^ 


int  ^ 3 A'  A'  B'  B 
2tt 


(2.  23) 


R 


(c ) The  strength  of  the  interaction 

We  calculate  the  interaction  energy  between  two  charges, 
at  r^  and  at  r^  , on  Al,  treating  the  substrate  as  a jellium, 

with  r = 2.  If  we  neglect  screening,  we  obtain  AE^  by  taking  the 


int 


external  potential  to  be 


Acp°(r) 


r - r 

/V 


r - r 

V A.5' 


(2.  251 


To  take  screening  into  account,  we  note  that  since  the  interaction 
energy  involves  interaction  between  electron -hole  pairs  of  wave  num- 
ber q 2 k^  , the  effect  of  screening  can  be  approximated  by  using  a 


static  surface  dielectric  function  e„  at  q - 2 k , so 

o F 


that 


AE. 


AE. 


int 


int 


(2.  2b) 


For  r^  - 2,  the  bulk  dielectric  function  e corresp>onding  to 
q = ^ value  of  1.  15,  ^ ^ only  slightly  different  from  1. 


1 1 


The  surface  dielectric  function  e is  not  known  at  this  time  but  it  is 
clear  that  it  will  also  differ  only  by  very  little  from  unity.  In  what 
follows,  we  set  - 1. 

The  interaction  energy  (from  eq.  (Z.  23),  (2.  14),  (2.9)  and 

(2.  25))  is  plotted  in  Fig.  II.  For  definiteness,  we  have  taken  r^  = 2, 

z , - z „ - (2  k_)  ^ (=  0.  52  atomic  units)  and  Z,  = Z„  = 1.  For  other 
ABF  12 

values  of  and  z^,  the  amplitude  of  the  interaction  is  obtained  by 

scaling  by  the  ratio  of  the  product  I(z  ) and  I(z  ),  which  are  given 
in  Table  1. 

III.  Interaction  Energy:  Tight- Binding  Model 

In  this  section,  we  deal  with  a tight-binding,  simple  cubic 
semi -infinite  solid  with  a (100)  face.  The  adatoms  A and  B,  each 
with  a single  energy  level  E , sit  on  top  of  the  substrate  atoms  at 
surface  sites  0 and  n respectively.  The  Hamiltonian  of  the  system 
is  taken  to  be 

H = Hq  + Hj  + XH^  . (3.  1) 

Hq  is  the  Hamiltonian  for  the  solid,  and  has  nearest  neighbor  matrix 
elements  only,  which  are  denoted  by  T,  so  that  the  bandwidth  of  the 
solid  is  12  T,  centered  at  zero  energy.  H^^  is  the  Hamiltonian  for 
the  adatoms,  with  eigenvalue  E . XH,  represents  the  interaction  be- 
tween  the  adatoms  and  the  solid.  Its  only  non-zero  matrix  elements 


I 


4 


12 


are  V = V - V = V = V.  ),  is  treated  as  a small  parameter, 
AO  Bn  OA  nB  ^ ’ 

A/  ^ /V 

but  will  eventually  be  taken  as  unity.  Screening  is  neglected  since, 

as  we  saw  in  Sec.  II,  its  effect  is  small. 

To  obtain  the  indirect  interaction  energy  between  adatoms  A 

and  B,  we  have  to  perform  a perturbation  calculation  to  fourth -order 

in  X . More  conveniently,  we  can  use  the  method  of  canonical 
(13) 

transformation. 


~ -iXS  iXS 
H = e He 


S is  chosen  so  as  to  remove  to  lowest  order  in  X the  matrix  elements 

A/ 

of  which  transfer  electrons  between  the  adatom  and  the  substrate 

atoms : 


H^  +i[H^  + H^,  S]  = 0 


(f  alH_lf'a) 
(f  a|s|f'a'>  = -i— ^ 


for  i ^ £ 


and  <fa|s|fa')  = 0,  ^ 

where  If)  and  |f  are  eigenstate  s of  Hq  > which  are  Slater  deter- 
minants of  single -particle  eigenstates  Ik)  . |a)  and  \a')  are 

eigenstates  of  H and  are  of  the  form  I n , n ) where  n , n _ = 0 
1 A B A B 


With  this  < hoice  of  S,  we  have 


(3.  5) 


/N/  2 ^ 

H = Hq  + + X . 


where  the  relevant  matrix  element  of  is  now  given  by 

|H.|f  a>  =y  S (f'a|H  |f"a,">  (f"a"|H  |l  a> 


^ f"a 


(3.  6) 


We  have  two  cases,  depending  on  whether  the  unperturbed 

eigenvalue  E_  is  greater  than  or  less  than  the  Fermi  level  E . For 
u F 

the  case  E^  < E^,  , eq.  (3.  6)  describes  a virtual  excitation  from  the 

ground  state  |f  , a)  to  the  excited  state  "a")  , in  which  the  electron 

on  one  of  the  adatoms  is  transferred  to  the  conduction  band  above  the 

Fermi  level.  The  state  a)  is  formed  when  another  electron  in 

the  conduction  band  below  the  Fermi  level  is  then  transferred  to  the 

adatom,  creating  an  electron-hole  pair.  The  situation  is  similar  for 

the  case  E.  > E^  . 

0 F 

The  single  particle  eigenstates  |k)  of  are  of  the  form 


sin(k  (m  + a))  m ) , 
z z ~ 


where  |m)  is  a localized  wavefunction  centered  at  m . m is  the 

2 

component  of  m perpendicular  to  the  surface,  and  m =0  for  sur- 
f\j  z 

face  sites.  N is  the  number  of  lattice  atoms,  a is  the  lattice  con- 


stant and  k - (p  k ). 

~ jo  z 


14 


By  second-order  perturbation  theory,  treating  the  term  \ 
as  the  perturbation,  we  have 


AE. 


int 


i6|v^“  y 
2 


N 


2 i(P  -P 

sin  (k'  a)  sin  (k"  a ) e ~ 

k',k"  ^ ^ 


• na 


k F 


1 ...  ' 


^k'-S  V-^0  ""k'-V 

/V  A/ 


1 


(3.8) 


When  = E , eq.  (3.  8)  becomes  singular  and  we  have  a 
resonance  . This  corresponds  to  virtual  transfer  of  electrons  be- 
tween the  adatom  and  the  conduction  band,  leading  to  a very  strong 

interaction.  If  E„  is  not  too  close  to  E , and  E is  near  the  bot- 
U r r 

tom  or  the  top  of  the  energy  band,  so  that  the  Fermi  surface  is  nearly 
spherical,  eq.  (3.  8)  can  be  evaluated  in  the  same  way  as  (2.  16)  to 
give  (putting  X = 1 ), 


Re 

(2tt)^ 


-+ 


1 


^0  \"-^o 


f 2 i(p^ - p" ) • na 
k'^k"  e - - 


z z 


Ta^  [2k  (k'  -k")  + k^  -k"^l 

L F X X z z J 


-L Jvil_ 

4 2 2 

T (E^-E^) 


k a V ^ cos(2,k  _na) 


X 2T 


(3.  9) 
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sects  the  zone  boundary  and  has  saddle  points.  The  method  used  in 


Sec.  II  to  evaluate  the  interaction  energy  would  not  be  valid,  because 

(14) 


of  the  different  geometry. 


IV.  Interaction  Involving  Surface  Bands 

If  the  substrate  has  a surface  band  that  is  partly  filled,  virtual 
excitations  and  de -excitations  of  electrons  between  empty  and  occupied 
states  in  the  band  due  to  the  presence  of  the  adatoms  will  contribute  to 
the  interaction  energy. 

We  assume  that  the  Fermi  surface  for  the  surface  band  is 
cylindrical,  so  that  the  energy  denominator  in  (Z.  3)  is  given  by 


^0-^i 


1 ,2  /2. 

— T (P  - P ) 


2 m 


(q  - 2q  • p)  , 


(4.  1) 


2 m 


where  m ' is  the  effective  mass,  p and  p'  are  two-dimensional 
wave  vectors  parallel  to  the  surface  and 


q = p - p 

/V 


(4.  2) 
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The  wavefunctions  'I' (r)  are  periodic  in  the  x-  and  y-direc- 
tions  but  damped  in  the  z-direction  on  both  sides  of  the  surface. 
Analogous  to  (Z.  11),  the  matrix  element  ( i jcp  | 0 ) is  given  by 


(i  jcc  |0  > = e^S  ■ S V(q)  , 


(4.  3) 


where  V(q)  has  the  same  form  as  (2.  11)  except  that  the  wavefunctions 

/N/ 

f(z;  k ) are  replaced  by  damped  waves.  As  before,  the  only  relevant 
z 

matrix  elements  in  the  summation  in  (2.  3)  involve  momentum  transfer 

q =»  2k  ^ where  k^  is  the  Fermi  wave  number  for  the  surface  band, 

so  that  V(q)  can  be  replaced  by  V , evaluated  at  q = 2k^  . Equa- 

OF 

tion  (2.  3)  then  becomes 


^ (0|cpg|i)(i|cpg|0) 

^ V.  I~F 


E - E. 
o 1 


+ c,  c. 


2-*  |Vq|" 


k,  q, spin 


iq  • R 
e ^ /v 

1 2 

-y  (q  - 2q  • p) 


* 1 I 2 

2m  |VJ 


»)‘  /■ 


2 iq  • R 
d q e G(q)  , 


(4.  4) 


whe  re 


G(q)  = 


■2(q  -2qp, 


(4.  5) 


A is  the  surface  area  of  the  substrate  and  the  factor  2 takes  spin  into 


account.  As  in  Sec.  11,  we  differentiate  G(q)  with  respect  to  q and 
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evaluate  the  leading  tei-m  in  (q  - Zk^).  We  find  that  the  most  singular 

r 

term  is  given  by 


so  that 


-^G(q)  = Iq  - 9(q  - 2kp  . 


^ ■ r>  ->  cos(2k'  R) 

j2  lq•R^,^  -,2  F 

d q e ~ G(q)  = -2n  


(4.7) 


m'lVgl^A^  cos(2  k'  R) 


(4.  8) 


-2  -5 

The  interaction  falls  off  like  R , as  compared  to  the  R 
behavior  for  the  bulk  band.  Thus,  when  a surface  band  is  present  and 
the  Fermi  level  lies  inside  the  band,  its  contribution  to  the  interaction 
will  dominate  at  large  distances. 

To  determine  the  magnitude  of  the  interaction,  we  need  the 
value  of  Vq  . We  assume  that  the  wavefunction  'l'(r)  is  periodic  in  x 
and  y and  damps  exponentially  away  from  the  surface.  Since  typi- 
cally the  surface  wavefunction  damps  much  more  quickly  on  the  vacuum 
side,  we  just  cut  it  off  on  this  side  for  simplicity.  The  normalized 
wavefunction  is  approximated  by 


ip  • R Kz  „ , 
e X,  ~ e 9 (-z)  , 


(4.  9) 
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where  K is  a positive  constant.  For  a bare  charge  at  r = (0,  z ) 

~ o 

we  have 


2K  / 

A J 


V(q)  = ^ 5e^^R^+  (z  - z 

o 

z S 0 


0 


/ / J„< 


2K  / 2Kz  / 2 ? -1/2 

2n  / dze^""^  / J _(qR)[R  + (z-z^)*^]  RdR 


0 


/ 


4ttK  / 2Kz  -1 


e q e 


dz 


-qz 


4nK  e 


Aq  2K  + q 


At  q - 2 k , we  have 
F 


V , ILK  e 
0 


- 2 k_  z 
F o 


A k^(k^,  +K) 


(4,  10) 


(4.  11) 


If  kj,  is  sufficiently  large,  we  may  neglect  screening  as  in  Sec.  11, 
so  that  the  interaction  energy  is  given  by 


AE. 


m"  cos(2k^R) 


int 


k'^^  (M(k^/K))^ 


(4.  12) 


R 


This  is  shown  in  Fig.  IV,  with  r =2  and  z = (2k  )"^  as  in  Fie  11 

8 o F • 

For  simplicity,  we  have  taken  m*  = 1,  k'  = k and  (k'  /K)  = 1. 

AX*  X* 
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V.  Discussion 


Grimley  has  derived  expressions  for  the  indirect  interaction 
between  chemisorbed  atoms  on  bcc  and  fee  crystal  surfaces.  Assuming 
that  the  Fermi  surface  is  spherical,  he  obtained  two  cases,  depending 
on  the  direction  of  the  vector  R joining  the  surface  sites  to  which  the 
adatoms  are  bonded.  In  the  critical  direction,  i.  e.  R • t =0, 

/>✓  f^Q 

where  is  one  of  the  three  primitive  translation  vectors  which  is 
an  element  of  symmetry  lost  in  going  from  the  infinite  to  the  semi- 
infinite crystal,  he  obtained  an  oscillatory  interaction  that  varies  as 
-5 

R . In  other  directions,  he  had  again  an  oscillatory  interaction,  but 
- 3 

varying  as  R . We  have  rederived  his  formulas  for  the  interaction 

-5 

energy  but  found  that  it  varies  as  R for  all  directions  of  R , in 
agreement  with  the  result  obtained  in  Sec.  II. 

Einstein  and  Schrieffer calculated  the  indirect  inter- 


action between  adatoms  on  "cubiurn”  for  various  values  of  the  parame- 
ters V,  E and  E . They  found  the  interaction  to  be  oscillatory, 
and  over  their  range  of  calculation  they  fitted  it  by  an  exponential 

form.  They  also  conjectured  that  asymptotically  the  interaction  amp- 

-5 

litude  should  vary  as  R , but  felt  that  the  asymptotic  power  behavdor 

was  not  applicable  to  the  relatively  small  separations  they  considered. 

-5 

CXir  (jwn  theory  does  indeed  yield  an  R behavior  and  we 
find  this  asymptotic  form  fits  the  numerical  results  of  Einstein  and 
Schrieffer  rather  well.  Some  of  their  results  are  shown  in  Fig.  111. 


J 


The  agreement  with  the  present  work  is  very  good  in  one  case,  and 
somewhat  worse  in  the  other.  The  reason  is  presumably  that  in  the 
second  case,  the  Fermi  surface  is  quite  anisotropic,  and  replacing 
it  by  a spherical  surface  is  not  a very  good  approximation. 

Tsong,  using  field  ion  microscopy,  observed  that  the  inter- 
action between  Re  atoms  on  the  (110)  face  of  W is  oscillatory,  with 
a period  of  about  5 A,  and  falling  off  slowly  with  distance,  approxi- 
mately as  R \ The  interaction  can  be  fitted  by  an  expression  of  the 
form 


AE.  ^ = -C 
int 


cos(Ak  R) 
o 


with  C = 0.  95  X 10  and  Ak  = 0.  665  atomic  units.  Because  of  the 

o 

complexity  of  the  W Fermi  surface,  we  have  not  attempted  a detailed 
calculation  for  this  system.  We  wish,  however,  to  make  the  follow- 
ing observations  which  we  believe  to  be  relevant  for  this  system. 

The  W Fermi  surface  consists  of  an  electron  "jack"  at  the 
center  of  the  Brillouin  zone  and  hole  "octahedra"  at  the  corners  H, 

with  portions  of  rather  flat,  parallel  faces  normal  to  the  (111)  and 

0-1  (15-17) 

equivalent  di  rections,  approximately  Ak  - 1.0  A apart. 

For  virtual  excitations  of  the  substrate  electrons  of  wave  vector  q 
in  the  (111)  direction  (parallel  to  the  substrate  (110)  face),  these  flat 
faces  provide  large  regions  around  the  Fermi  surface  for  which  the 


energy  denominator  in  (2.  3)  is  nearly  zero,  and  hence  give  the  major 


contribution  to  the  interaction  energy. 


We  perform  a model  calc:ulation  with  a Fermi  surface  defined 

by  two  exactly  parallel  planes,  Ak  apart,  perpendicular  to  the  (111) 

direction  (taken  as  the  x-direction).  The  zone  boundaries  in  the 

other  two  directions  are  adjusted  so  as  to  give  the  correct  volume  in 

k-space.  This  gives  a rectangular  box  of  dimensions  Ak  X k^  X k^  , 
4-1  0-1 

with  Ak  = 1.  0 A and  k^  0.  69  A . The  energy  denominator  in 

' 3.  8)  is  given  by  (k  ^ - k ' ^)/2m  where  m*  is  the  effective  mass. 

XX 

After  integrating  over  k and  k'  , (3.8)  becomes 

z z 


AE  . 


4lv|4a^ 


int 


(2n)^  (E„-E  )^ 

F o 


2 m k / dk  dk  dk  dk  ^ 


; -11  \ 

X y y y q(q  - 2k^) 


Ivl  a^m  ''k^  / 2 sin  ~1 


2tt  (E„-E  ) 
F o 


f -1-ln 

q + Ak  1 

/ 2q 

q - Ak  1 

e dq 


I i4  6 * , 2 

|V  a m k 


2tt  (E  -E  ) 
F o 


JT cos( Akx) 

Ak  1^1 


(5.  2) 


The  pair  correlation  function  for  the  adatoms  is  given  by 

g(R)  = exp[- AE.^^(R)/kgT]  . (5.3) 

At  high  temperature  (k  T » AE.  ,(R))  the  spatial  dependence  of  the 

o 1 nt 

interaction  energy  can  be  obtained  by  averaging  (5.  2)  over  all  angles. 

2-1 

Since  ^E^^^  is  proportional  to  (xy  ) , we  expect  that  the  angular 
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average  should  vary  as  R 


At  low  temperature,  the  adatoms  are  "frozen"  at  sites  of  mini- 
mum energy,  corresponding  to  y = 0 and  cos(Akx)  = 1,  Thus, 


Ivl"*  ^ ^ 

' ‘ O COS{Akx) 

2tt^(E„  - E Ak  h ! 

F o 


cos  ( Akx) 


(5,4) 


This  has  the  same  form  as  (5.  1),  falling  off  as  lx|  and  with  a 
period  of  oscillation  given  by  (2rr/Ak)  - 6.  3 A,  compared  to  the 
observed  period  of  5 A . 


The  value  of  C'  in  (5.4)  can  be  estimated  as  follows.  The 
substrate,  which  is  simple  cubic,  is  assumed  to  have  the  same 
density  as  W.  Since  W is  bcc,  with  lattice  constant  = 3.  16  A , 
this  gives 


° - 1 / 3 4 

a = 3.  16  A X 2 = 2.  51  A . 


(5.  5) 


=:=  (17)  0-1 

m Rs  1 ' and  k « 0.  69  A . The  value  of  (E_  - E ) is  deter- 
o F o 

mined  by  assuming  that  the  energy  band  of  W is  formed  by  the 

broadening  of  the  5d  electrons,  with  a bandwidth  = 1 2 T == 

(17  ) 

10.  5 eV,  centered  at  the  5d  atomic  level.  Since  there  are  four 

5d  electrons  in  W,  E„  is  at  0.  35  T below  the  center  of  the  band. 

r 

E is  assumed  to  be  at  the  5d  level  of  atomic  Re.  The  5d  levels  of 
o 


W and  Re  are  at  -0.  684  and  -0.781  Rydberg  respectively. 


(18) 

These  give  ( E_  - E ) = 0.  7 1 eV  (see  Fig.  V).  To  obtain  a value  of 
“ F o 

_ 2 

C'  0.  95  X i 0 a.  u.  <-s  in  (5.  1 ),  we  need  V - 1.  67  eV  = i.  8 1 . 

I'  iS  IS  quite  reasonable,  in  view  of  the  fact  that  Einstein  and 
Schrieffer  estimated  that  V/T  be  roughly  in  the  range  of  3 to  5 for 
most  chemisorption  systems. 

in  Tseng's  experiment,  the  interaction  energy  is  of  the  order 
of  30  meV  (~  350°K)  and  the  experiment  was  performed  at  320°K,  so 
that  the  temperature  is  in  the  intermediate  range.  Neither  of  the 
extreme  cases  considered  above  is  therefore  applicable,  but  the 
experimental  observations  appear  to  be  qualitatively  consistent  with 
the  effects  discussed  in  this  paper. 
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FIGURE  CAPTIONS 
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Fig.  I 


Fig.  II 


Effective  one-electron  potential  (jellium  model; 

r = Z). 
s 

Interaction  energy  between  two  bare  charges  near  the  jellium 

surface  (r  = Z).  The  charges  are  at  a distance  (Zk_)  ^ from 
s F 

the  jellium  edge. 


Fig.  Ill  Comparison  of  the  interaction  energy  in  the  tight-binding 
model  with  the  calculations  of  Einstein  and  Schrieffer 
(open  circles). 

Fig.  IV  Interaction  energy  involving  surface  bands. 

Fig.  V Schematic  diagram  for  the  energies  in  the  Re/W  system. 

All  energies  are  in  eV. 
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